Topological entanglement entropy relations for multi phase systems with interfaces 
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We study the change in topological entanglement entropy that occurs when a two-dimensional 
system in a topologicaUy ordered phase undergoes a transition to another such phase due to the 
formation of a Bose condensate. We also consider the topological entanglement entropy of systems 
with domains in different topological phases, and of phase boundaries between these domains. We 
calculate the topological entropy of these interfaces and derive two fundamental relations between 
the interface topological entropy and the bulk topological entropies on both sides of the interface. 
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In gapped two dimensional systems, the ground state 
entanglement entropy Se of a spatial domain with 
boundary of length L scales with L as Se — oiL — 7 + • • • . 
Here a is a constant and the dots represent contributions 
that vanish in the limit L — 00. The additive constant 
—7 is called the topological entanglement entropy [HH]. 
It is independent of the detailed geometry of the spa- 
tial region (depending only on its topology) and its value 
is robust against perturbations of the system's Hamilto- 
nian, changing only in phase transitions. In fact, one may 
express 7 explicitly in terms of the quantum dimensions 
da of the various types of topological excitations that 
may occur in the system For a region with the topology 
of a disk, 7 = \ogD, with D = \/X)a '^a- Topologi- 
caUy trivial phases have just one topological sector with 
quantum dimension equal to 1 and hence have D = 1 
and 7 = 0, while all topologicaUy nontrivial phases have 
7 > 0. Hence 7 has become a popular indicator of the 
presence of topological order. In fact the size of 7 is an 
indicator for the complexity of the topological order or 
topological symmetry present in a system. In particular, 
D and the number of topological sectors together deter- 
mine whether a system hosts non-Abelian anyons, since 
such anyons have quantum dimensions greater than 1. 

It is clearly of interest to know how 7 changes in phase 
transitions In Refs. [3H5], we have developed a method for 
the description of topological symmetry breaking transi- 
tions between topological phases, caused by condensation 
of a bosonic excitation. In such transitions, the value of 
D always decreases. Our first aim here is to calculate 
the corresponding increase in bulk topological entropy in 
terms of a number q characterizing the transition. After 
that, we obtain a formula for the topological entropy of a 
spatial interface between two phases related by conden- 
sation. Finally we consider multilayer systems. 

Let A and U denote topological phases of the same sys- 
tem, where U is obtained from A by topological symme- 
try breaking through a Bose condensate. We will abuse 



notation and call the quantum groups, unitary braided 
tensor categories or fusion algebras of these phases A and 
U as well, as long as no confusion can be expected. 

Since 7^ = log Da and = legist/, the topological 
entropy increases by 

-7c/ = log^r- (1) 

We will now derive an expression for this number in terms 
of data directly related to the condensation transition. 

The features of the description of topological symmetry 
breaking are as follows (see Refs. [S^5]). We assume that 
a bosonic excitation of phase A condenses. The new, con- 
densed, topological phase has particle like (localized) ex- 
citations whose fusion and braiding is described by quan- 
tum group U . However, the condensed phase also allows 
for confined excitations which pull physical strings in the 
condensate. In a system with domains of both phases, the 
confined excitations are expelled to the boundary of the 
condensed domain. On this boundary there is a well de- 
fined set of fusion rules for both confined and unconfined 
excitations, described by a fusion algebra T. While T al- 
lows for consistent (associative) fusion, it will not usually 
allow for a consistent implementation of the braid group, 
exactly because the confined excitations are attached to 
the boundary by strings. Our scheme therefore involves 
three levels, or two steps. We start with the unbroken 
phase described by the quantum group A, which by the 
condensate gets broken to an algebra T. A subset of the 
T-sectors corresponds to unconfined excitations. These 
are described by a quantum group U 

The breaking step is described by a branching of topo- 
logical charge sectors of the A-phase into sectors of the 
intermediate fusion algebra T, 

a— >^n*t (restriction). (2) 

t 

Here, the are integer coefficients. We can think of this 
decomposition of A-sectors as resulting from restriction 
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of representations of the quantum group A to a subalge- 
bra which is not broken by the phase transition. Con- 
versely the lift of the fields from T to A is defined by 



(lift), 



(3) 



where n° — ria- Crucially, restriction must commute with 
fusion in the A and T theories, and preserve the quantum 
dimensions of the sectors, so da — J2t ''^a^t- 

For any topological sector t G T, we can define a num- 
ber q related to the transition from A to U hy 



dt 



(4) 



where the di which appear are quantum dimensions in 
the A and T theories. In the appendix, we show that q 
does not depend on the choice of t. In fact, we prove that 
q — D\/D'^. In addition, we will soon argue on physical 
grounds that — DaDjj (cf. Eq. ([o])), so that we have 

q = Da/Du. (5) 

These facts about q may also be proved rigorously with 
the use of advanced mathematics, either in the framework 
of modular tensor categories or using von Neumann al- 
gebras [7] , but we will only use elementary arguments and 
physical intuition here. Given Eq. ([5]), we may conclude 
that the change in topological entanglement entropy be- 
tween the phases is simply given by 



- 7c/ = log q- 



(6) 



so q characterizes the phase transition in the sense that it 
determines the accompanying loss of topological entropy. 

In Ref. [S], we argued that, when A and U are topologi- 
cal phases which can be described by conformal field the- 
ories (CFTs), quantum group symmetry breaking from 
A to U can occur if there exists a conformal embedding 
A <Z U oi the chiral algebras U and A of those CFTs. 
Therefore, we call q the quantum embedding index [T^ . 

A simple example of the breaking scheme discussed 
above occurs for A = SU{2)4, giving U = SU{3)i upon 
condensation of the spin 2, or weight 4 particle in the 
SU{2)4 theory. This case correspond to the conformal 
embedding of the SU {2)4 affine Lie algebra into SU{3)i 
and it is extensively discussed in Ref. [5j. In table |lj we 
give the sectors of the two theories with their quantum 
dimensions, as well as those of the intermediate algebra 
T. One may readily verify that indeed: 
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2. (7) 



While q is an integer in this case, usually this will not be 
true. For example when A = S'[/(2)io and U = SO{^)i_ 
(see Ref. [5] for details), one finds that q = ?> + ^/2>. 
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TABLE I: Topological sectors A with their quantum dimensions 
d/^ and topological spins /ia , as they feature in the breaking of an 
SU{2)4 theory after condensation of A = 4 excitations. Decom- 
position of A-sectors into T-sectors is indicated. The spins of U 
particles are the same as for the corresponding yl-particles. Spins 
of confined T-particles are ill-defined. 



Let US now consider a geometry as indicated in Fig. [T] 
where we have taken a sphere with the left hand side in 
an unbroken A phase and the right hand side in the con- 
densed U phase, separated by the interface described by 
T. We have also divided the sphere into two regions / 




FIG. 1: 

(shaded gray) and // for which we wish to calculate the 
topological entanglement entropy. Since the entangle- 
ment entropy of complementary subsystems is equal, we 
have Sej = Seji- As the interface geometry is the same 
for both regions, the topological entanglement entropies 
must also be equal. Region I consists of two disks, one in 
the A-phase domain and one in the [/-phase domain, so 
the total topological entropy here equals —7^ — ■^ij. We 
also would like to calculate the entropy of the annular 
region II directly. The topological entropy for an annu- 
lus was calculated by Kitaev and Preskill in Ref. [T] and 
was found to equal twice the entropy for a disk. How- 
ever, their calculation was done for a system in a single 
topological phase and does not directly apply here as the 
region // contains both phases. Still, if we narrow the 
annulus // down to a circle so that only the states on the 
interface contribute, we would expect to obtain topolog- 
ical entropy — 7t- These arguments suggest a relation 
between the bulk and interface topological entropies. 



7A + 7c/ = 277 



which implies that 



DaDu = D^, or Da/Dt = Dt/Di 



(8) 
(9) 
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We observed this relation already Ref. [S], but now we 
have given it a physical interpretation. From Eq. ^ we 
also see that Da/Dt = Dt/Du = ^g, so D is reduced 
by a factor of ^ in each of the two steps of the topolog- 
ical symmetry breaking process. 

Now consider a two-layer geometry of two concentric 
surfaces as indicated in the left picture of Fig. [2] with the 
inner and outer layers in phases Ai and A2, resp. Assume 
that, in some domain, a Bose condensate forms in a chan- 
nel that couples the Ai and A2 phases' degrees of free- 
dom. We then have a situation as indicated in the right 




FIG. 2: Side view of double layered system with A ^ A^® A2 
where on the right a condensate has formed (in gray area). 
The situation corresponds to the coset construction if A\ = 
G,A2^H and U = G/H. 

picture of Fig. [2] which leads to the entropy relation: 
lA^+lA^+lu = 27T, or Da.Da^Du = D^. (10) 

The geometries depicted in Fig. [2] can in fact represent 
the coset construction of CFTs [5j, if we take Ai — 
G,A2 = H and U = G/H. In that case we have 

iG+lH+lG/H ^ "^St or DgDhDc/h = D^, (H) 

underscoring the role of T in the coset construction. 

Next imagine that we take away the disc labeled A2 
from the uncondensed area, yielding the situation of 
Fig. [3] We introduced this setup in Ref. W to describe 




FIG. 3: Single layered system bordering a two layer conden- 
sate (gray region). In this situation the central charges of the 
two domains can be different and gapless degrees of freedom 
can then exist on the domain boundary. 

transitions and interfaces between phases at different cen- 
tral charge c, such as those occurring between different 
fractional quantum Hall conductance plateaus. In Bose 
condensation transitions, c is conserved, even though D 
changes. Adding or removing a layer however also adds 
or subtracts the central charge of that layer. One now 



obtains a modified entropy relation; clearly 7/ is lowered 
by ^-nd therefore the same has to hold for 777. Since 
T has not changed, we find that the interface topolog- 
ical entropy is now —277^ + 7^2- We can heuristically 
interpret this increase of interface topological entropy as 
follows. Because of the difference in the central charges 
on either side, the interface now allows for excitations de- 
scribed by a chiral CFT with the same value of c as the 
broken layer. In the simplest cases, this is just an A2 type 
CFT (More detail will appear in [in])- The topological 
entanglement entropy —7 normally receives a (negative) 
contribution from every topological sector of the theory, 
in gapped bulk phases, where the topological excitations 
are localized. However the excitations described by the 
A2 CFT are gapless and delocalized and it appears they 
do not count toward 777. 

As an example of a two level system consider A — 
Ai®A2 ~ I sing (g) I sing and assume a condensate forms 
in the bosonic (e, e) channel. Table [TT| lists the fields of the 
intermediate T algebra with their quantum dimensions, 
and the lifts of the fields into A. The fields (cr, 1) and 
(1, cr) are confined and the remaining fields have Z2 (8)^2 
fusion rules. We have Da = 4, Dt = 2\/2 and D{j = 
2, and q — Da/Djj — 2. Removing one of the Ising 
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TABLE II: The breaking scheme of the A = I sing ® I sing theory, 
after formation of the (e, e) condensate: A => T = Z2 ^ Ising => 
U = Z2 22. The total quantum dimensions equal Dj\ = 4, = 
2V2 and Du = 2. 

layers, as in Fig. |3j we obtain a boundary between a 
theory of Ising type and one of the toric code universality 
class. Such a situation can be realized microscopically in 
the Kitaev hexagonal models (cf. TT]). The interface 
between regions will carry a CFT of Ising type and the 
value of 777 will then be reduced by \og{Dising) = log(2). 
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Appendix: Relations for q and D 

We will now derive that the quantum embedding index 
q as defined in Eq. ([d]) is independent of t and satisfies 
q — D\/D^. We describe both the A and T theories at 
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the level of fusion algebras. For A, we write the fusion 
rules as a X b = N^^^c and define the fusion matrices Na 
by {Na)l — iV^f,. For T we use the same notations, but 
with indices r, s,t, . . . rather than a, b,c, . . .. The fusion 
coefficients for A and T satisfy iV^^ = Nj^^, and for A we 
also have N^^, = N^^. 

We write the restriction as a map B : A T, given 
by Ba = J2a "a^- ^^^^ given by the transpose map 
B^ : T ^ A. We require that fusion and restriction 
commute, so Ba x Bb = B{a x b), ov more explicitly, 

^7V>* -^n^iV;,. (12) 

c r.s 

We also have conservation of quantum dimension under 
the restriction, J^t "■a'^t = '^a, as well as = ng. 

In both A and T, the quantum dimensions satisfy the 
fusion rules, e.g. for A, we have 

dad, = N^,d,. (13) 

As a result, the vector uja defined by uja = ^a^a-O- is 
a common Perron- Frobenius eigenvector of all the Na, 
satisfying the eigenvalue equations NaiOA = daUiA- Sim- 
ilarly, wt — J2t^tt is a common eigenvector of the iVj, 
despite the fact that these do not necessarily commute. 

We need the behavior of the Perron-Frobenius eigen- 
vectors uiT and UJA under the restriction B and the lift 
B'f. We have that, B^ut = UJA, since, in components, 
this is just the conservation of quantum dimension under 
the restriction. We will now show that also 

BujA — q^T (14) 

For this, let us define the map M : T ^ T hy 

Af^^A/„ = ^<iVt (15) 

a at 

M is the fusion matrix of the restriction of J^a ^- "^^^ 
matrix elements of M are all strictly positive, since for 
any r,s € T there is some t G T with N^,. > and 
also any t d T occurs in the decomposition of some 
a £ A. Therefore the eigenspace of M corresponding to 
the Perron-Frobenius eigenvalue da is non-degenerate 
and spanned by ujt- Hence, if we can show that 

MBujA = iJ2'ia)BuJA (16) 

a 

then it follows that Buja = Awt for some A e C. In 
components, Eq. ( [16^ reads 

a,c,s,r a,b 

To prove this, one can massage the left hand side using 
iV^j = and ria — n^. Then one apphes (|l2| , uses 
that iV|g = N^f^ = iV^^ and finally one applies ( |13[ ). Some 
intermediate results are given below, 

J2 dcnlnlK, = dcnlKNl, = ^ d^A 

a,c,s,r a,c,s,r a,h,c 

= 'YdcNli,nl = "Ydadbn\- 

a,b,c a, 5 



Hence Buja = ^ojt- Writing this in components, we find 



A= * ° =q 

dt 

with q as defined in Q . This shows that q is independent 
of the choice of t. Since q is the quotient of the quantum 
dimension of a field i in T (or mU <ZT) by the quantum 
dimension of its lift in A, we call it the quantum embed- 
ding index. Another reason for this terminology is the 
similarity to the embedding index introduced by Dynkin 
in his work on the classification of subalgebras of com- 
pact Lie algebras. This index satisfies a similar formula, 
where the quantum dimensions should be replaced by the 
indices of the Lie algebra representations. 

From Buja = <-^t and B^ujt — qui a, it follows that 
B^B and BB^ satisfy 



B^ Buja = Q^A and BB^ ujt = quJT (18) 



From this we obtain that q — D\/D^ as follows. 



D\ = {a*\a*) = {t*\B^B\t*) ^q{t*\t*) ^qD^. (19) 



Combining this relation with Q, one also gets ([s]). 
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